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Abstract 

The properties of ultrathin films have been studied within the framework of Ising model and 
the method of random-field interactions. It is shown that the Curie temperature is inversely 
proportional to the number of layers. Critical exponent v has been obtained and it is shown 
that it does not depend on the type of crystalline lattice. 

As a characteristic, that allows to evaluate the smallness of magnetic objects, one can take the 
length of the spin-spin correlation It is used to describe the spin fluctuations near the critical 
temperature. In the case when the film thickness is less than the correlation length of bulk samples, 
there appear various size effects. In particular, the decrease of film thickness can lead to a reduction 
of phase transition temperature: according to measurements on films of iron, nickel and cobalt [l}|8], 
Tc decreases with decreasing of film thickness A'^ (here N is the number of monolayers). The exper- 
imentally observed dependence of the Curie temperature on the thickness Tc{N) can be described 
by the relation introduced in [6jj7|: 

def 

where Tc {N — t- oo) = T^o - temperature of bulk sample, cq and N - constants. The argument A is 
related to the critical exponent of the spin-spin correlation u as follows: A = l/u [o]. And, as was 
shown in 10-12 , the constants and argument of the dependency ([T]) essentially depend on the type 



and symmetry of the crystal lattice of sprayed film and its substrate. 

In this paper we assess the influence of film thickness and its crystal structure on the temperature 
of magnetic phase transition. To solve this problem we use the following model. 



1 Model 

• Ultrathin film is composed of N infinite monolayers; 

• the interaction fields h between spin magnetic moments of the atoms are distributed randomly, 
and the interaction is realized only between the nearest neighbors; 

• spin magnetic moments are oriented along an axis Oz (approximation of the Ising model) and 
are equal to mo in magnitude. 

According to [ 13] the distribution function for random interaction fields /i on a particle located 
at the origin (in n-th monolayer) is defined as: 

Wn{h) = J ^ ^-"^ Kj (rj, nij) j ]^ Fn (mj) 5 (r^ - rjQ)drjdm.j, (2) 

where 6{x — xq) - Dirac delta function, /i„j = hnj (nij, rj) - field produced by atoms with magnetic 
moments nij, located at the nodes with coordinates Vj in n-th monolayer. Fn (nij) is the distribution 
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function for the magnetic moments, which in the approximation of the Ising model for a ferromagnet 
can be represented as follow: 

Fn (m,) = {anjS {Oj) + /3„,(5 - ii)) 6 (mo,) . (3) 

Here 9j is the angle between rrij and Oz-axis, a„j and /3„j - relative probabilities of the spin 
orientation along {9j = 0) and against {9j = vr) Oz-axis respectively; moj - magnetic moment of the 
j-th magnetic atom. 

Now we introduce symmetric notations a_|_„ a„ and a-n 1 — an and consider the set kj 
of nearest neighbors of an arbitrary atom numbered as j. Let denote z = dim kj as its coordination 
number. Then let Q^j be the set of all possible products of a±n (where n £ kj), which contains z 
various factors with different values of n (therefore dimQfc^^ =2^). The elements of $7^^, are marked 

as u>f UJi Ilnefc ~ Y\'{u=i- i^ek ) '^ilv fact, £ is a number of binomial permutation of a±n 
that forms $7. If L is such a set of binomial permutations, on which Q was build, then we can define 

def 

a similar set M of elements Mi: Mg = XlnGfc = ^0 Snefc ^l^"^" ~ -'-I- 

In the approximation of nearest neighbors and the direct exchange interaction between magnetic 
atoms, the equation ^ can be represented as: 

where {kj) is a sample of n atoms of the total number of z nearest neighbors of j-th atom, Ji^ are 
the constants of exchange interaction (that may differ between different monolayers or even inside 
the same monolayer between different sorts of atoms). 

General equation that determines the average relative magnetic moment in n-th monolayer 

is 

^in = J tanh (j^^ Wn{H)dH. (5) 

Replacing in expressions for and Mi all a±n on their average values {a±n) = (1 =b fin)/'^ and 
substituting (Q into (|5]), one can obtain the equations that determine /i„ in each monolayer: 

/=o 

Z^fc=0 '-'21,2 2^ \ t 

H-n — 1^1=0 ^Zn,n 2^"." 2^k = '"'Zn-l.n 2"^"-l." (6) 

^Z7l,n+l (l+/Jn+l)''(l-A'n+l)"'"'"+^~'' 4-„„v, / 2 (;+A:+r) - (z:^- 1 +2:n,n +2:„,„+i ) ^ 



Mat — Z^i=0 ^ZN,N 2^'v,iV X 



where z„^„ is the number of nearest neighbors in the n-th layer, Zn-i^n is the number of nearest 
neighbors of the atom in (n — l)-th layer, located in the n-th layer; i„„ = J„„m„/Jiimi, in-i,n = 
Jn-i,ninn-i/Jiimi, in,n+i = Jn,n+imn+i/ Jiimi, t = k T/Jnmi. Using one can study the 
dependence of the average magnetic moment of the film on its temperature and thickness. 



2 The Curie temperature of ultrathin films 

Fig. jljshows the temperature dependence of the average relative magnetic moment (m) = {/J,j) /N 
films of different thickness. From the graphs it follows that the decrease in the number of mono- 
layers leads to a reduction in the average number of nearest neighbors and, consequently, to lower 
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FCC(lll) 




Figure 1: The dependence of the average magnetic moment (m) on the temperature for different 
crystallographic planes of the FCC lattice in the films of different thickness. On the inset on the left 
it is shown an increase in the average number of nearest neighbors [z) in films of different crystalline 
structures, depending on the thickness of the A^. 



the transition temperature Tc- Moreover, the temperature dependence of {m) = {m{T)), as well 
as the position of the Curie point, determined not only by the type of crystal lattice, but also by 
the crystallographic orientation of the plane of the film growth. The above-mentioned feature is 
related to the difference in the number of nearest neighbors between atoms in the same monolayer, 
and atoms located in adjacent layers. For example, in a material with a FCC lattice grown on 
(lOO)-plane, each atom in the monolayer has 4 neighbors, and 4 more neighbors in the monolayer, 
which is located nearby. When the film of the same material grows on (lll)-plane, the number of 
nearest neighbors changes to 6 and 3, respectively. 

The dependence of the Curie temperature Tc on the thickness of films with different crystal 
structures, that have been calculated by solving (|6]), is shown on the Fig. |2] From the illustration it 
is clear that Tc essentially depends on the type of crystal lattice and is almost independent of the 
crystallographic orientation of the surface of the film, except for the area of small thickness (2 — 6 
monolayers). This effect can be explained by the fact that in the area of small amount of monolayers 
the mean number of neighbors differs notably. 

The calculation of relative change of the phase transition temperature with the growth of film 
thickness (in logarithmic scale) is shown in Fig. [s] The dependence e = £{N) can be approximated 
by expression ([T|. The results can be compared with the experimental data (see the table). Within 
the measurement errors and the accuracy of approximation, the calculated values of A and the 
critical exponent of the spin-spin correlation u are close to the experimental that obtained on films of 
Ni/Cu{lll) and Ni/Cu{100). We note that the argument A does not depend on the crystallographic 
orientation of the surface (Ano ~ Am), nor the type of lattice (A for FCC and BCC lattices differs by 
0.9%). In addition, the critical exponent of the spin-spin correlation {u) = j = 0.68, calculated in 
this framework, is close to the value 1^= 0.63, obtained by using renormalization group calculations 



in three-dimensional Ising model 14 15 



3 Conclusion 

The model of randomly interacting atomic magnetic moments at its relative simplicity provides two 
valuable results. Firstly, it allows to assess the influence of the thickness of ultrathin film on the 
temperature of magnetic phase transition, and, secondly, to establish a power-dependence of the 
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Figure 2: The dependence of the reduced Curie temperature Tq on the thickness of ultrathin film 
(A^ is the number of monolayers) for different crystal lattices (FCC and BCC) and different crystal- 
lographic orientations of the surface. Compare with the inset in Fig. [T] 




Figure 3: The logarithmic dependence of relative change of the Curie temperature £{N) on the film 
thickness N (in monolayers) for FCC (on the left) and BCC lattice. 



Table of theoretical and experimental values of argument A and the critical exponent i/: 



Film 


I (monolayers) 


A 




Link 


Ni/Cu(lll) 


1 - 8 


1.48 ±0.20 


0.68±0.09 




10 1 


Ni/Cu(lll) 


1 - 10 


1.44 ±0.20 


0.70±0.10 




11 j 


Calculation: FCC (111) 


2-10 


1.43 


0.70 




Ni/Cu(100) 


4-26 


1.42 


0.70 




12] 


Calculation: FCC (100) 


2- 26 


1.48 


0.68 
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relative change of the Curie temperature of N: e{N) ~ N^^'^, where u is the critical exponent of 
the spin-spin correlation. The obtained regularities and numerical values of critical exponents agree 
well with experimental data (6{[7{ |10}jT2] as well as with RG calculations in three-dimensional Ising 



model 14 15]. The magnetic properties of the lattice (in particular, the independence of the choice 



of crystallographic orientation of the surfaces and the type of lattice) are consistent with the general 
concepts of critical scaling. 
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